Refining the transcriptional landscapes for distinct clades of virulent phages infecting Pseudomonas aeruginosa.

The introduction of high-throughput sequencing has resulted in a surge of available bacteriophage genomes, unveiling their tremendous genomic diversity. However, our current understanding of the complex transcriptional mechanisms that dictate their gene expression during infection is limited to a handful of model phages. Here, we applied ONT-cappable-seq to reveal the transcriptional architecture of six different clades of virulent phages infecting Pseudomonas aeruginosa. This long-read microbial transcriptomics approach is tailored to globally map transcription start and termination sites, transcription units, and putative RNA-based regulators on dense phage genomes. Specifically, the full-length transcriptomes of LUZ19, LUZ24, 14-1, YuA, PAK_P3, and giant phage phiKZ during early, middle, and late infection were collectively charted. Beyond pinpointing traditional promoter and terminator elements and transcription units, these transcriptional profiles provide insights in transcriptional attenuation and splicing events and allow straightforward validation of Group I intron activity. In addition, ONT-cappable-seq data can guide genome-wide discovery of novel regulatory element candidates, including noncoding RNAs and riboswitches. This work substantially expands the number of annotated phage-encoded transcriptional elements identified to date, shedding light on the intricate and diverse gene expression regulation mechanisms in Pseudomonas phages, which can ultimately be sourced as tools for biotechnological applications in phage and bacterial engineering.

Phage proteins target and co-opt host ribosomes immediately upon infection.

Bacteriophages must seize control of the host gene expression machinery to replicate. To bypass bacterial anti-phage defence systems, this host takeover occurs immediately upon infection. A general understanding of phage mechanisms for immediate targeting of host transcription and translation processes is lacking. Here we introduce an integrative high-throughput approach to uncover phage-encoded proteins that target the gene expression machinery of Pseudomonas aeruginosa immediately upon infection with the jumbo phage ΦKZ. By integrating biochemical, genetic and structural analyses, we identify an abundant and conserved phage factor ΦKZ014 that targets the large ribosomal subunit by binding the 5S ribosomal RNA, and rapidly promotes replication in several clinical isolates. ΦKZ014 is among the earliest ΦKZ proteins expressed after infection and remains bound to ribosomes during the entire translation cycle. Our study provides a strategy to decipher molecular components of phage-mediated host takeover and argues that phage genomes represent an untapped discovery space for proteins that modulate the host gene expression machinery.

Exploring the transcriptional landscape of phage-host interactions using novel high-throughput approaches.

In the last decade, powerful high-throughput sequencing approaches have emerged to analyse microbial transcriptomes at a global scale. However, to date, applications of these approaches to microbial viruses such as phages remain scarce. Tailoring these techniques to virus-infected bacteria promises to obtain a detailed picture of the underexplored RNA biology and molecular processes during infection. In addition, transcriptome study of stress and perturbations induced by phages in their infected bacterial hosts is likely to reveal new fundamental mechanisms of bacterial metabolism and gene regulation. Here, we provide references and blueprints to implement emerging transcriptomic approaches towards addressing transcriptome architecture, RNA-RNA and RNA-protein interactions, RNA modifications, structures and heterogeneity of transcription profiles in infected cells that will provide guides for future directions in phage-centric therapeutic applications and microbial synthetic biology.

SAPPHIRE.CNN: Implementation of dRNA-seq-driven, species-specific promoter prediction using convolutional neural networks.

Data availability is a consistent bottleneck for the development of bacterial species-specific promoter prediction software. In this work we leverage genome-wide promoter datasets generated with dRNA-seq in the Gram-negative bacteria Pseudomonas aeruginosa and Salmonella enterica for promoter prediction. Convolutional neural networks are presented as an optimal architecture for model training and are further modified and tailored for promoter prediction. The resulting predictors reach high binary accuracies (95% and 94.9%) on test sets and outperform each other when predicting promoters in their associated species. SAPPHIRE.CNN is available online and can also be downloaded to run locally. Our results indicate a dependency of binary promoter classification on an organism's GC content and a decreased performance of our classifiers on genera they were not trained for, further supporting the need for dedicated, species-specific promoter classification tools.

Transcriptional Organization of the Salmonella Typhimurium Phage P22 pid ORFan Locus.

Many phage genes lack sequence similarity to any other open reading frame (ORF) in current databases. These enigmatic ORFan genes can have a tremendous impact on phage propagation and host interactions but often remain experimentally unexplored. We previously revealed a novel interaction between phage P22 and its Salmonella Typhimurium host, instigated by the ORFan gene pid (for phage P22 encoded instigator of dgo expression) and resulting in derepression of the host dgoRKAT operon. The pid gene is highly expressed in phage carrier cells that harbor a polarly located P22 episome that segregates asymmetrically among daughter cells. Here, we discovered that the pid locus is fitted with a weak promoter, has an exceptionally long 5' untranslated region that is instructive for a secondary pid mRNA species, and has a 3' Rho-independent termination loop that is responsible for stability of the pid transcript.

A Grad-seq View of RNA and Protein Complexes in Pseudomonas aeruginosa under Standard and Bacteriophage Predation Conditions.

The Gram-negative rod-shaped bacterium Pseudomonas aeruginosa is not only a major cause of nosocomial infections but also serves as a model species of bacterial RNA biology. While its transcriptome architecture and posttranscriptional regulation through the RNA-binding proteins Hfq, RsmA, and RsmN have been studied in detail, global information about stable RNA-protein complexes in this human pathogen is currently lacking. Here, we implement gradient profiling by sequencing (Grad-seq) in exponentially growing P. aeruginosa cells to comprehensively predict RNA and protein complexes, based on glycerol gradient sedimentation profiles of >73% of all transcripts and ∼40% of all proteins. As to benchmarking, our global profiles readily reported complexes of stable RNAs of P. aeruginosa, including 6S RNA with RNA polymerase and associated product RNAs (pRNAs). We observe specific clusters of noncoding RNAs, which correlate with Hfq and RsmA/N, and provide a first hint that P. aeruginosa expresses a ProQ-like FinO domain-containing RNA-binding protein. To understand how biological stress may perturb cellular RNA/protein complexes, we performed Grad-seq after infection by the bacteriophage ΦKZ. This model phage, which has a well-defined transcription profile during host takeover, displayed efficient translational utilization of phage mRNAs and tRNAs, as evident from their increased cosedimentation with ribosomal subunits. Additionally, Grad-seq experimentally determines previously overlooked phage-encoded noncoding RNAs. Taken together, the Pseudomonas protein and RNA complex data provided here will pave the way to a better understanding of RNA-protein interactions during viral predation of the bacterial cell.IMPORTANCE Stable complexes by cellular proteins and RNA molecules lie at the heart of gene regulation and physiology in any bacterium of interest. It is therefore crucial to globally determine these complexes in order to identify and characterize new molecular players and regulation mechanisms. Pseudomonads harbor some of the largest genomes known in bacteria, encoding ∼5,500 different proteins. Here, we provide a first glimpse on which proteins and cellular transcripts form stable complexes in the human pathogen Pseudomonas aeruginosa We additionally performed this analysis with bacteria subjected to the important and frequently encountered biological stress of a bacteriophage infection. We identified several molecules with established roles in a variety of cellular pathways, which were affected by the phage and can now be explored for their role during phage infection. Most importantly, we observed strong colocalization of phage transcripts and host ribosomes, indicating the existence of specialized translation mechanisms during phage infection. All data are publicly available in an interactive and easy to use browser.

Introducing differential RNA-seq mapping to track the early infection phase for Pseudomonas phage ɸKZ.

As part of the ongoing renaissance of phage biology, more phage genomes are becoming available through DNA sequencing. However, our understanding of the transcriptome architecture that allows these genomes to be expressed during host infection is generally poor. Transcription start sites (TSSs) and operons have been mapped for very few phages, and an annotated global RNA map of a phage - alone or together with its infected host - is not available at all. Here, we applied differential RNA-seq (dRNA-seq) to study the early, host takeover phase of infection by assessing the transcriptome structure of Pseudomonas aeruginosa jumbo phage ɸKZ, a model phage for viral genetics and structural research. This map substantially expands the number of early expressed viral genes, defining TSSs that are active ten minutes after ɸKZ infection. Simultaneously, we record gene expression changes in the host transcriptome during this critical metabolism conversion. In addition to previously reported upregulation of genes associated with amino acid metabolism, we observe strong activation of genes with functions in biofilm formation (cdrAB) and iron storage (bfrB), as well as an activation of the antitoxin ParD. Conversely, ɸKZ infection rapidly down-regulates complexes IV and V of oxidative phosphorylation (atpCDGHF and cyoABCDE). Taken together, our data provide new insights into the transcriptional organization and infection process of the giant bacteriophage ɸKZ and adds a framework for the genome-wide transcriptomic analysis of phage-host interactions.